GENERAL TRANSFORMATIONS BETWEEN THE HEUN 
AND GAUSS HYPERGEOMETRIC FUNCTIONS 



GALINA FILIPUK 



Faculty of Mathematics, Informatics, and Mechanics, University of Warsaw 
Banacha 2, 02-097 Warsaw, Poland 
E-mail: filipuk@mimuw. edu.pl 



RAIMUNDAS VIDUNAS 



Faculty of Mathematics, Kobe University 
Rokko-dai 1-1, Nada-ku, 657-8501 Kobe, Japan 
E-mail: vidunas@math.kobe-u.ac.jp 



Abstract. It is known that Heun's functions (or differential equations) can be reduced to Gauss 
hypergeometric functions by rational changes of its independent variable only if its parameters, 
including the fourth singular point location parameter and the accessory parameter, take special 
values. We present all hypergeometric-to-Heun transformations with two or three free continuous 
parameters up to fractional-linear transformations. 

1. Introduction. The ordinary differential equations for the general Gauss hypergeometric 
function and the Heun function are second order Fuchsian differential equations on the Rie- 
mann sphere with 3 or 4 regular singular points, respectively. Any second order Fuchsian 
equation with 3 or 4 singularities can be transformed to them by Mobius transformations 



1.1. The hypergeometric equation. The differential equation for the general Gauss hyper- 
geometric series 
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has the form 
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It has regular singularities at z = 0, 1 and oo. The local exponents are given by 

0, 1 — catz = 0; 0, c — a — b at z — I; and a, 6 at z = oo. 

Recall that if A is a local (characteristic) exponent at a regular singularity z = t of the 
Fuchsian differential equation, then in a neighborhood oi z — t there is a local solution 
y{z) = {z — t)'*' J2T=o ^k{z — t)^ . Accordingly, the hypergeometric equation ([2]) with general 
a, 6, c has the following local bases of solutions: 



at z = : 
at z = 1 : 

at z = oo : 



2F1 
2F1 



a, b 
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z^-' 2F1 



1 + a — c, 1 
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c — a— 6 
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-b ib,l + b-c 
1 — a + 
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c — a, c — 6 
1 + c — a — & 
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z 



1 - z 



The local exponent differences at singular points are (up to a sign) equal tol — c, c — a — 6 
and a — 5, respectively. In this paper we assume that the local exponent differences are 
non-negative. The information about singularities and local exponents is encoded in the 
Riemann scheme for ([2]): 

1 00 

pI a z 

1 — c c — a — b b 

Recall that Fuchsian equations with 3 singularities are defined uniquely by their singularities 
and local exponents. This is not generally true for classes of Fuchsian equations with more 
singularities; accessory parameters are then needed to specify a Fuchsian equation. 

The following Pfaff and Euler fractional-linear transformations 2, Th. 2.2.5] can be 
applied to 6 local solutions: 



2F1 



a, b 



= (l-z)-%Fi 



= (I-z)-SFi 



a, c — 
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c~ a, b 
c 



z 



z - 1 
z 



(1-z)' 



: — a — b 



z - 1 
c — a, c — b 
c 



(3) 



In total we have 6 x 4 = 24 different hypergeometric solutions for a general hypergeometric 
equation; they are referred to as 24 Kummer's solutions. For integer values of the parameters 
a,b,c or the local exponents, the structure of 24 solutions may degenerate [M]. The 24 
solutions are related by a group of transformations that permutes the singular points or 
interchanges the local exponents at the singular points. The group is a semidirect product 
of the symmetric group S3 and (Z/2Z)'^, but the interchange of local exponents at infinity 
does not change the hypergeometric function. These transformations of the hypergeometric 
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equation can be presented as transformations of Riemann's P-synibols; for example 

P< a z\ ^ {1 - zY'^'-^'P I c-b z 

— cc — a — 6 b ) — ca + 6 — cc — a 

r 1 oo 

= P\ a I- z 

\ c — a — b 1 — c b 

For a basic theory of hypergeometric functions see [5], [H Vol. I]. 

1.2. The Heun equation. The Hewn equation is defined by 

cPy{x) f c d a + b — c — d + 1\ dy{x) abx — q 



dx^ ' \x ' x — 1 ' x~t J dx ' x{x — l){x — t)^^ ^ ^ 

The equation has 4 regular singular points: x — 0, x — 1, x = t, x — oo. The parame- 
ters a, b, c, d determine the local exponents at them, while the parameter q is an accessory 
parameter. The Riemann scheme of the Heun equation is 

r 1 t oo 

P < ax 

— c 1 — d c + d — a — b b 

A local solution at x = is defined by a power series which coefficients satisfy a three-term 
recurrence relation. For comparison, the 2F1 local solutions of the hypergeometric equation 
satisfy a similar two-term recurrence relation. In the following we shall denote the solution 
of (dl with the local exponent and the value 1 at x = by 



Hn 



(5) 



t a, b 
q c] d 

This function reduces to the Gauss hypergeometric function ifd = a + b'-c+l and 
q abt; notice the reduction of equation ([4]) to (l2|) then. The function is identical to the 
constant 1 if afe = and q — 0. Note that the parameters a, b are symmetric and give the 
local exponents at x = c», whereas the parameters c and d are not interchangeable and 
determine the non-zero local exponents at x = and x = 1, respectively. 

Similarly as with the hypergeometric equation, permutation of the 4 singular points and 
interchange of local exponents at them act on Heun's equation by changing the parameters, 
including the location t of the fourth singular point and the accessory parameter q. The 
acting group is a semi-direct product of S4 and (Z/2Z)'*. The group has 4! x 2^ = 384 
elements, but the interchange of the local exponents at x = 00 does not change Heun's 
equation. Correspondingly, there are 384/2 — 192 solutions expressible in terms of Heun's 
function ((S]); they are presented in |8j. Here we present a basic generic structure of the 192 
solutions. 

Permutation of the singular points x=l,x = t, x = ooor interchange of local exponents 
at these points, are realized by two-term fractional-linear transformations of Heun's func- 
tion. First of all, we have transformations that do not move the singular points, but only 
interchange the local exponents at x = 1 and x — t. These transformations neither change 
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the argument x nor the location parameter t: 



Hn 



a, h 
c; d 



Hn 



q-cid~l)t 



c-\-d—a — b 



Hn 



91 



a-d+1, 5-ff+l 
c; 2 - d 
c + (i — a, c + d — b 
c: d 



(6) 



(1-x) 



1 



c+rf— a— b 



Hn 



c - a + l, c-b+l 
c\2~d 



where qi — q ~ c{a + b — c — d), q2 =q~c{a + b — c — d + dt — t). Permutations of the 
three singular points x — 1, x^t, x = oo are realized by the following fractional- linear 
transformations : 



Hn 



a, b 
c: d 



fl/t 



Hn 

(1-x 



a, b 

c; a + b — c — d+l 



' \{act-q)/{t-l) 



1 - - 
t 



Hn 



1/(1 -t) 
(q-ac)/(t-l) 



1 - - 
t 



Hn 

(1 - j;)""Hn 



1 - t 
ac — q 
1 - 1/t 
ac — q/t 



a, a — d + 1 
c; a — 6 + 1 
a, c + fi — 6 
c; a — 6 + 1 
(1 - i)a; 



S - 1 
X 

X t 



(7) 



a, c + d — 6 

c; d a; — t 

a, a — d + 1 
c; a + 6 — c — d+l 



(i- l)a; 



i (a; - 1), 

Combined with the transformations of 4 functions above, we get 6 x 4 = 24 expressions for 
the local solution ^ at x = 0. Like with the hypergeometric equation, we generally have 
two local solutions at each of the 4 singular points expressible in terms of Heun's function. 
Here are general bases of local solutions at a; = 0, .t = 1, a; = f and x = oo: 
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Hn( 

a;"^Hn( 



qi = (c- l)(a + 6-c-(i + (it + l), 

(72 — ab — q — {d — l){a + b — ct — d + 1), 

q^ = ab — q/t + (c/t — c — d){a + b — c — d), 

qi = q/t + a{a-b/t- c - d + d/t+l), 

95 = q/t + b{b~ a/t- c- d + d/t+1). 



l-x 



(8) 
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Each of these 8 functions can be expressed in 24 ways by fractional-hnear transformations 
©-(III), giving the 192 solutions in 8 . Any 3 of the 8 functions are related by a linear 
three-term connection formula (since the order of Heun's equation is 2), though their coeffi- 
cients are not known in general (unlike Rummer's solutions of the hypergeometric equation). 
The group of transformations of Heun's equation can be represented by two-term transfor- 
mations of Riemann's P-symbols, or by the action on the 4 signed local exponent differences 
(1 — c, I — d, c + d — a — b — a) . A transformation can multiply some of the local exponent 
differences by —1, and permute them. 

The Heun equation contains a large number of interesting special cases, in particular, 
the Lame equation, which is of considerable importance in mathematical physics j4l llOj . 
The Heun equation appears in problems such as diffusion, wave propagation, magneto- 
hydrodynamics, heat and mass transfer, particle physics and cosmology of the very early 
universe. Heun functions arc much less well understood than hypergeometric functions. In 
particular, no general elementary integral representation of Heun functions is known. 

2. Pull-back transformations of differential equations. For an ordinary linear differ- 
ential equation a rational pull-back transformation has the form 

z ^ ^{x), y{z) ^ Y{x) = e{x) y{^{x)), (9) 

where ip{x) is a rational function and 9{x) is a radical function, i.e., a product of powers 
of rational functions. The term 9{x) is called a pref actor in the terminology of Maier in [7]. 
Geometrically, transformation ([9]) pull-backs a differential equation on the projective line 
Fl to a differential equation on the projective line P;^, with respect to the finite covering 
(/? : PiJ, — > determined by the rational function (p{x). Here we let F]., P], denote the 
complex projective lines with rational parameters x, z respectively. In [1], [15] these pull- 
back transformations are called RS -transformations. 

In this paper we are interested in general pull-back transformations of the hypergeometric 
equation to the Heun equation. These transformations typically give two-term transforma- 
tion formulas 

between the Heun and Gauss hypergeometric functions. These identities give reductions of 
Heun equations (with special values of the parameters) to the better understood hypergeo- 
metric functions. It is generally desirable to have expressions of Heun's and especially Lame 
functions in terms of more elementary functions. 

We classify all such transformations with 2 or 3 free continuous parameters, up to 
fractional-linear transformations ([6|)-(l7|) of hypergeometric and Heun's functions. A 
similar problem is considered by Maier in [7j, where all transformations ()10|) without the 
prefactor 9{x) and coming from pull-back transformations of their differential equations are 
classified, with any number of free parameters. In this paper, we allow prefactors 9(x). There 
is an overlap between our and Maier's lists. We plan to classify all transformations (I10|l with 
1 or free continuous parameters in future papers. The remainder of this section gives a 
technical background of our project. 

The relevant coverings ip : F^. ^ F\ typically branch only above the 3 singular points 



Hn 



if{x) (10) 
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z — 0,z— l,z — ooof the hypergeometric equation. We describe the behavior of singularities 
under general pull-back transformations in Subsection 12.11 below. The restriction that the 
transformed differential equation must have only 4 singular points turns out to be very 
limiting. Pull-back ransformations ([9]) between hypergeometric equations or between Heun 
equations are classified similarly. In [1] , [15j the same methodology is applied to obtain pull- 
back transformations to isomonodromic 2x2 Fuchsian systems corresponding to algebraic 
solutions of the Painleve VI equation. 

2.1. General action of a pull-back transformation. The following definition is introduced in 
[12] . An irrelevant singularity for an ordinary differential equation is a regular singularity 
which is not logarithmic, and where the local exponent difference is equal to 1. An irrelevant 
singularity can be turned into a non-singular point after a suitable pull-back transforma- 
tion ^ with ip{x) = X. For comparison, an apparent singularity is a regular singularity 
which is not logarithmic, and where the local exponents are integers. Recall that a loga- 
rithmic point is a singular point x — t where there is only one local solution of the form 
{x — t)^(l -\- X]fc°=i Oikix — t)'')- For us, a relevant singularity is a singular point which is not 
an irrelevant singularity. 

A pull-back transformation of a Fuchsian equation gives a Fuchsian equation again, usu- 
ally with more singularities. We are interested in pull-back transformations of the hyperge- 
ometric or Heun's equation to Heun's (or hypergeometric) equation. Since any second order 
Fuchsian equation with 3 or 4 singularities can be transformed to the hypergeometric or 
Heun's equation by a Mobius transformation, the essential problem is to look for pull-back 
transformations to a Fuchsian equation with 4 (or 3) singular points. Relevant singular 
points and local exponent differences for the pull-backed equation are not influenced by the 
prefactor 9(x). The number of relevant singularities is determined by the covering ip alone. 
Finding suitable coverings is our essential problem. 

The behavior of regular singular points of a Fuchsian equation under a pull-back trans- 
formation is summarized in the following results from 12J. 

Lemma 2.1. Let </? : P;J, ^ be a finite covering. Denote a Fuchsian equation on by Hi, 
and the pull-back transformation of Hi under (0) by H2. Let S denote a point on V^., and 
let k denote the branching index of tp at S . 

1. The local exponents for H2 at S are equal to kai -\- 7, ka2 + 7, where: 

• ai, a2 are the local exponents for Hi at ^(S) G P^; 

• J is the local exponent of the radical function 6{x) at S . 

2. Lf the point ip{S) is non- singular for Hi, then the point S is not a relevant singularity 
for H2 if and only if k — \ (i.e., the covering does not branch at S). 

3. Lf the point ip^S) is a singular point for Hi , then the point S is not a relevant singularity 
for H2 if and only if the following two conditions hold: 

• The point ip{S) is not logarithmic. 

• The local exponent difference at ip{S) is equal to 1/fc. 

Proof. The first statement is mentioned in the proof of [T^ Lemma 2.4]. The other state- 
ments are parts 2 and 3 of [121 Lemma 2.4]. ■ 
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Lemma 2.2. Let Lp ^ ¥\ be a finite covering of degree d. Denote a set of 3 points on 
by IS.. If all branching points of ip lie above A, then there are exactly d + 2 distinct points 
on P;J, above A. Otherwise there are more than d + 2 distinct points above A. 

Proof. Part 1 of [121 Lemma 2.5]; follows from the Hurwitz formula 5, Corollary IV. 2. 4]. ■ 

Suppose we start with a hypergeometric equation Hi on P^; let A denote the set {0, 1, oo} 
of its singularities. It follows from the above lemmas that to keep the number of singularities 
of the pull-backed equation low, we should allow the branching points of (f only above A 
(because the outside branching points would be revelant singularities, and we would get more 
than d + 2 distinct points above A), and we should restrict some local exponent differences 
of Hi to the values 1/fc, where k is an integer. Recall that a covering of P"'^ by a Riemann 
surface is a Belyi covering [TT] if it branches only above a set of three points. It follows that 
relevant coverings Lp will typically be Belyi coverings. Non-Belyi coverings can occur only in 
the following rather degenerate situations. 

Lemma 2.3. Suppose we have a pull-back transformation (0) of a hypergeometric equation 
Hi to Hennas equation, and there is a branching point o/ (yS : P;J, — > Pj that does not lie above 
{0, 1, oo} G ¥\. Then one of the following statements holds: 

• Hi has actually only 2 relevant singular points; 

• Two of the 3 local exponent differences of Hi are equal to 1/2; 

• Hi has a basis of algebraic solutions. 

Proof. Let d denote the degree of ip. By Lemma [221 there are at least d + 3 distinct points 
above {0, l,oo} S P^- At most 3 of them can be singularities of the pull-backed equation 
(because the outside branching points are singularities). Therefore we have at least d non- 
singular points above {0, 1, oo}. If we restrict a local exponent difference of Hi to 1/k with a 
positive integer fc, we have at most d/k non-singular points above that point. We can choose 
k — 1, but then the point is non- logarithmic only if it is not a relevant singularity by [121 
Lemma 5.1]. Otherwise we can restrict two local exponent differences to 1/2, or restrict all 3 
local exponent differences to 1/k, l/£, 1/m with l/k + l/£ + l/m > 1 and k,£,m integers. If 
l/k-\-l/£-\-l/m = 1 and we optimally arrange the non-singular points above {0, 1, oo}, we 
have only those d non-singular points above {0, l,oo}, contradicting Lemma 12.21 The case 
1/k + l/£ + 1/m > 1 gives algebraic solutions. ■ 

Note that the monodromy group of Hi in the three exceptional cases of Lemma 12.31 is. 
respectively, a cyclic group, a dihedral group, or a finite group. 

Apart from these special cases, the rational argument (p{x) of a transformation ()10p 
between the Heun and hypergeometric functions will define a Belyi covering, and the pull- 
backed equation will have 4 singular points and d — 2 non-singular points above the singular 
locus {0, 1, oo} of the hypergeometric equation, where d is the degree of (p{x). For compar- 
ison, typical transformations between hypergeometric functions have 3 singular and d — 1 
non-singular above the same locus. 

We are interested in pull-back transformations to the Heun equations with at least two 
continuous parameters. Therefore we want to restrict none or just 1 local exponent difference 
of a starting hypergeometric equation. We may ignore the possibilities with k — 1, because 
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non-logarithmic hypergeometric equations with a local exponent difference 1 have only one 
continuous parameter by [121 Lemma 5.1]. 

Referring to a Belyi covering ip of degree d, we denote the branching pattern above the 
three critical points by three partitions of d separated by the equality sign, just as in [12] . 
For example, a degree 3 Belyi covering can have the branching pattern 2 + 1 = 3 = 2 + 1, 
meaning one fiber with a single point with the maximal branching index, and two fibers 
with a simple branching point and a non-branching point. The total number of elements in 
the three partitions must be equal to d + 2. 

To describe transformation of local exponent differences under a pull-back transforma- 
tion, we extend the notation (ai,/3i,7i) (a2,/32,72) in [H] used for transformations 
between hypergeometric equations. In particular, we write (ai,/?i,7i) (02, /?2, 72, <52) 
to denote a pull-back transformation of degree d of a hypergeometric equation with the 
local exponent differences ai , /3i , 71 to Heun's equation with the local exponent differences 
0^2, /32, 72, 1^2. The arrow follows the direction of the covering Lp -.V]. ^T\. The order of local 
exponent differences on both sides of the arrow is irrelevant for us, as we do not assign the 
local exponent to particular singularities with the notation. 

2.2. Adjusting the prefactor and other technical issues. Once we have a suitable covering 
(f -.Fl. ^ Pi and choose local exponents of the starting Fuchsian equation, the prefactor 
9{x) for a target pull-back transformation ([9|) is determined up to a constant multiple, and 
it is straightforward to derive it. For identification of hypergeometric or Heun solutions of 
both equations, the prefactor is normalized by the condition 9{0) = 1. A direct pull-back 
transformation with 9{x) — 1 usually gives a differential equation with no local exponent 
at the points above z = 00 equal to 0, whereas standard hypergeometric or Heun equations 
have local exponents at all points except x = 00. Besides, all regular points on P;^ must 
have the local exponents and 1, hence the local exponents at irrelevant singularities must 
be shifted to 0, 1 rather than to 0, —1. Irrelevant singularities may occur above non-infinite 
z-points if the local exponents at them were restricted to 0, —1/k with k a positive integer. 

The prefactor must shift the local exponents at the points above z = 00 and all irrelevant 
singularities so that we would have the canonical local exponent value everywhere except 
possibly at a; = 00. The prefactor does not change the local exponent differences or the 
characteristic sum of all local exponents. (Recall that the sum of local exponents is equal to 
1 for the hypergeometric equation and to 2 for the Heun equation by the Fuchs relation). 
The prefactor has the form 9{x) — Yiii^ — Ci)^^*, where (t^ are all the points where the 
local exponents need to be shifted, and is the local exponent at Ui to be shifted to 0. 
The local exponents at x = 00 then shift by the sum of all ^^'s. It is convenient to use 
Riemann's P-notation to keep track of changes of local exponents. Intermediate P-symbols 
would display irrelevant singularities as well. 

The prefactor is not needed, (i.e., 9{x) = 1) if there is only one point a: = 00 above 
z = 00, and there are no irrelevant singularities above other z-points. The rational function 
ip{x) defining the covering is then a polynomial. 

Maier [7] classified all transformations between hypergeometric and Heun equations with- 
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out a prefactor. He obtained a list of five indecomposable transformations 

(a, /3, 7) ^ (a, a, 2/3, 27), (i, a, /?) ^ (i, a, 2a, 3/3) , (i, a, /3) ^ (a, a, a, 3^), 

V2' 3' / V2' 2' 3' V2'3' / V2'3'3' /' 

and two composite transformations (of degree 4 and 6): 

(i,a,/3) <^ (a,a,2/3) ^ (a, a, 2a, 4/3), ^ (i^Sa) ^ (1,1,1, 6a). 

Besides, more transformations without a prefactor are given for the degenerate Heun equa- 
tion with ab = q = 0. The function ip{x) does not to have to be a polynomial then, as the 
points above z — 00 immediately have a local exponent 0. 

To derive a two-term identity between hypergeometric or Heun solutions from a pull- 
back transformation of their Fuchsian equations, we choose a singular point on P], (above 
the singular locus on P^) of the pull-backed equation as a; = 0, and assign the label z = to 
its image on P^. The rational function (p{x) is determined once we label the other singular 
points on P], and Pj. The prefactor 6{x) is determined if additionally the local exponents for 
singularities on P^ are fixed. The formula identifies the local solutions a,t x — and z — 
with the local exponent and the value 1; it holds in a neighborhood x — once we set 

A choice of x = and local exponents at z = determines a two-term identity up to 
fractional-linear transformations ([3]), (I6|)-(l7|). If the local exponent difference at z = is 
restricted to 1/k with an integer k, the choices 0, 1/fc and 0, — 1/fc of local exponents give 
different identities. Changing the sign of the local exponent difference at z = basically 
gives an identity between the other two local solutions (with non-zero local exponents) at 
X — and z = 0. For transformations between hypergeometric functions, this situation 
is captured by ,12, Lemma 2.3]. Here is a reformulation for identities between Heun and 
hypergeometric functions. 

Lemma 2.4. Suppose that we have identity (|_?0p coming from a pull-back transformation 
between the corresponding hypergeometric and Heun equations. Then (p{x)^^^' ^ Kx^^'^ as 
x — > for some constant K , and the following identity holds: 



Hn 



1-fa-c, l-h6-c 
2 - c;d 



e(x)2F] 



l + A-C,l + B-C 
2-C 



where qi= q - [c - l)(a + 6 - c- d + dt + l) and Q{x) = 6{x) ip{xY'^ / Kx^^" . 

Proof. A straightforward identification of the other canonical local solutions of both equa- 
tions at a; = and z = 0. ■ 

If the local exponent difference at z = is an unrestricted parameter, changing its sign 
gives essentially the same two-term identity. More generally, the following choices of a: = 
give the same two-term identities up to the fractional-linear transformations and change of 
parameters: 

• the x-points with the same branching index and above the same point of P\ ; 

• the a;-points with the same branching index, if they are in different fibers with the 
same branching pattern, and either the local exponents at the corresponding z-points 
are the same, or the local exponent differences at both z-points are free parameters. 
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If the pull-backed equation is Heun's equation, its accessory parameter can be obtained 
at the latest stage, by considering power series expansions at a; = in a supposed two-term 
identity and comparing the first couple terms in the power series. 

Lemma 2.5. Suppose that we have identity (| iOp coming from a pull-back transformation 
between the corresponding hypergeometric and Heun equations, with 



as a; — > 0. Then q = ct \ n 



ABX" 



(x) = 1 + ^x + 0{x'^) 



C 



Proof. Expanding both sides of (|10p in the power series at x = gives 

1 + — a; + 0{x^) = 1 + ^ix + ^^\x + 0{x^). 
ct C 



In particular, if the covering ip{x) branches at a; = and the prefactor 0{x) is absent, 
then g = (because X — ^ = 0); see formulas (|16p and below. 

A pull-back transformation Q between hypergeometric and (or) Heun equations might 
fail to produce two-term identities between hypergeometric and (or) Heun solutions only if 
all singular points of the transformed equation lie above non-singular points of the starting 
equation. Lemma 12.31 then basically applies. So far we do not know an example of this kind. 
On the other hand, transformation identities between hypergeometric and Heun functions 
might formally exist without a pull-back transformation between their equations. For exam- 



ple, the linear function 2F1 ^ ^' ' 



z j can be formally transformed to any (hypergeometric or 
Heun) polynomial. Part 2 of [Ol Lemma 2.1] indicates that this situation can occur only if 
we start with a hypergeometric function actually satisfying a first order Fuchsian equation. 



2.3. Transformations between hypergeometric functions. Pull-back transformations between 
hypergeometric equations typically give algebraic transformations between hypergeometric 
functions of the form 



2F1 



a, h 



(x) 2F1 



A, B 
C 



(p{x] 



(11) 



The classical transformations were obtained by Gauss, Goursat, Riemann, Kummer. Here 
is an example of a cubic transformation, with one free parameter a: 



2F1I 



3a, i 
2a 4 



, / a, a + :t 



27a; 



(4a;- 1)3 



(12) 



The local exponents are transformed as (1/2, 1/3, a) < — (1/2, a, 2a), with a = 1/6 — 2a. 
Quadratic transformations have two free parameters. The transformation of local expo- 
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nents is (1/2, a, f3) 
2F1 

2F1 



(a, a, 2/3). Here are hypergeometric formulas: 




— 2F1 



1 



a, b 

+ b + 



4:X{1 - X) 




(13) 
(14) 

(15) 



The first two formulas are related by fractional- linear transformations ([3]), whereas (fTSll is 
not equivalent up to the fractional- linear transformations (on either P;J, or Pj), as noted by 
Askey [3] and Maier [9l Remark 4.1.2]. The dividing difference is the choice of the point 
X = 0: it is a non-branching point in (jl3p - (|14|) but a branching point in the last formula. 
But formula ifTS]) can be derived from (fT3|) by the following argument. The functions 



2a, a — b 
2a - 26 H 



1 



a, a -I- 



1 



(1 ~ 2xy 



^ 1 , (l-2a;)— '2F1I ' 
X ) \l + a — 

are among the 24 Kummer solutions of the differential equations for, respectively, the left- 
hand side and the right-hand side of (fT3|) . Therefore the two functions satisfy the same 
Fuchsian equation of order 2. We multiply both functions by x^°- , make the substitutions 
x I— > 1/x and 61-^a — 6-1-1/2 and obtain the two functions in ()15|) . up to a constant multiple 
on the right-hand side. Those two functions satisfy the same Fuchsian equation of order 
2, have the same value and the same local exponent at a regular singular point (with a 
non- integer local exponent difference in general), so they must be equal in a neighborhood 
of a; = 0, and p5|l follows. 

Pull-back transformations of the hypergeometric equations, and subsequently, algebraic 
transformations of the Gauss hypergeometric functions, are systematically studied and clas- 
sified in a series of papers by the second author; see |12j and references therein. Here is an 
example of a non- classical transformation (1/2, 1/3, 1/7) ^ (1/3, 1/7, 2/7): 

-1/28 



2F1 



( ^ 




19 




42 




42 






5 






7 







1 



19, 



343^2 , 
243 



16807 -j^ 
6561 I 




x^ (1 - x) (49a; - 81)^ 



4 (16807a;3 - 9261a;2 - 13851a; + 6561)3 I 

These identities can be verified by checking the power series at a; = 0. But the common 
region of convergence usually appears to be small. For example, the last identity does not 
hold at a; = 1 or a; = 81/49 for the standard analytic branches of 2F1 functions, as can be 
checked numerically. 

2.4. Quadratic transformations between the Heun and hypergeometric equations. By the 
Hurwitz formula, a quadratic covering (/s : P;]. —> has two ramification points. Application 
of a quadratic pull-back transformation to the general hypergeometric equation gives a Fuch- 
sian equation with 4 singular points, if the two ramification points lie above {0, 1, 00} C P^. 
Due to Mobius transformations on P;J,, we may assume that the points a; = 0, a; = l,a; = oo 
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are among those singular points, so that the puUbacked equation is Heun's equation. It 
foUows that quadratic transformations of general Gauss hypergeometric functions can be 
expressed in terms of Heun functions. Here are explicit formulas: 



Hn 



Hn 



-1 


2 
4a6 



2a, 26 
2c-l;a + 6-c+l 
2a, 26 

- 26 - 2c + 1 



c; 2a 
Hn 



1 

2 

2a6 



2a, 26 
c; c 



2F1 



— 2F1 



— 2F1 



a, 

c 
a, 6 

c 



x(2 



4x(l 



(16) 
(17) 
(18) 



They were first indicated by Kuiken in [6]. Other possible polynomials ip{x) for quadratic 
transformations between hypergeometric and Heun equations are l — x^, (1 — x)"^, {2x — 1)^. 
Fractional-linear transformations of P-symbols for the 192 Heun functions and the related 
Kummer's 24 hypergeometric functions give a set of another 30 rational functions of degree 
2 that transform the general hypergeometric equation to Heun's equations (with a prefactor, 
in general). The 30 rational functions are given in [6] in the context of the degenerate case 
a6 = g = 0. If we are not concerned about the prefactor, there is no distinction between 
quadratic transformations to degenerate or non-degenerate Heun equations. 

Like for hypergeometric quadratic transformations (|13p - (fT5|) . we have two different choi- 
ces for cc = 0: a branching point and a non-branching point. Accordingly, identities (|17p 
and (flS)) are related by fractional- linear transformations ([3]), ©-([T]), whereas identity (fT6| 
cannot be related to them by the fractional-linear transformations. To derive (|16p from (|17p , 
one can observe that the functions 



Hn 



-1 




2a 



2a, 26 
26 - 2c - 



1; c 



1 - x 



2F1 



a, b 
b-c 



1 



{i~xy 



satisfy the same Heun equation as both sides of (|17p . and have the same local exponent and 
value at a; = 1. Therefore they must be generally equal; formula (fTHl) is then obtained after 
the substitution 1 — x, cH^a + 6 — c+1. 

2.5. Heun-to-Heun transformations. Existence of quadratic and quartic transformations 
was pointed out by Erdelyi in [H Vol. 3] . Examples of these transformations are given by 
Maier in f9| Section 4]. Here are two alternative formulas of the quadratic transformations: 

1 



Hn 



Hn 



2a, 2a- 6 + 1 
6: 2a - 6 + 1 



2s-l 
2a6s-f 4qs(s- 



2a, 6 






6; 6 




(19) 



(20) 



where qi = (1 + s)'^q — 2abs. They are related by a series of fractional-linear transformations 
(and reparametrizations). The transformation of local exponents is given by the expression 
(1/2, 1/2, a, 13) < — (a, a, (3, (3); all choices for x — give two-term identities related by 
fractional-linear transformations (I6|)-(l7l). 

A quartic transformation can be obtained by composing two versions of the last quadratic 
transformation. In the composition, we have to restrict particularly s >— > l/2s, 6 1-^ 2a + 1/2 
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and a i— > 2a, h i-^ 2a + l/2 in the two versions. Remarkably, transformation of the parameters 



t and q simphfies greatly. After setting t 
transformation 

4a, 2a + i 
2a + i:2a + 



7(2.- 



Hn 



Aq 







ft 


a, a + i 






Hn 










2a+i; i 



1) in the composition, we recognize the 

Ux{x -l){x -t)\ 
(x2 - t)2 I 



as in [HI Theorem 4.2]. We note that a similar composition of quadratic transformations for 
Painleve VI functions is comparably simpler than the quadratic transformations as well [151 
Lemma 6.1]. The local exponents are transformed under the composition as 

(i, aj (i, i, a, a) (a, a, a, a). 

The composite degree 4 covering happens to be a Belyi covering. It is ramified over the 3 
points with the local exponent difference 1/2; the branching pattern is 2 + 2 = 2 + 2 = 2 + 2. 
Up to the Mobius transformation z 1/z, the starting Heun equation for the quartic 
transformation is a general Lame equation. 

Generally, Heun-to-Heun transformations of degree d must have at most 4 singular points 
and at least 2d — 2 non-singular points above the 4 singular points of the starting Heun 
equation, because the Hurwitz formula implies that a covering P], has at least 2d + 2 

points above any 4 points of . (An indicative situation is adding one fiber in the context 
of Lemma 12.21 ) It is not hard to see that the distribution of singular and non-singular points 
is impossible for d = 3 or d > 4, unless all local exponent differences of the starting equation 
are equal to 1/2 (or we allow it to have less singularities). We are then in a special case 
of the Lame equation; possible branching patterns indicate that the transformed equation 
has the same local exponents. Identification of solutions of the Lame equation with the 
local exponent differences (1/2,1/2,1/2,1/2) with the generic indefinite elliptic integral 
indicates that transformations of these Lame equations exist in any degree d; they arise from 
the isogenies between generic elliptic curves. Here is an example of a cubic transformation 
between respective Lame functions: 



Hn 



q{l~2sY 



Hn 



x{x + s{s-2)) 
{{l-2s)x + s^f 



(21) 



with s 7^ 0, 1/2, 1, 2. The branching pattern is a generic picture of 4 simple branching points 
in 4 different fibers. Note that there is no prefactor, even if the argument is not a polynomial. 



3. Transformations between the Heun and hypergeometric functions. Here we 
classify pull-back transformations ^ between hypergeometric and Heun equations with 2 
or 3 free continuous parameters (in the equations), up to fractional- linear transformations of 
their P-symbols. Subsequently, we present all corresponding Heun-to-hypergeometric trans- 
formation formulas (jlOp with 2 or 3 free continuous parameters, up to fractional-linear 
transformations ([3]), ©-([T]). 

Let Hi denote the starting hypergeometric equation, and d denote the degree of an 
assumed pull-back transformation ip : F], ^ Fl. 

If we assume 3 continuous parameters, we put no restriction on the local exponent 
differences of Hi. Therefore all points above z = 0, z — 1, z — oo will be singular for the 
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Local exponent differences 


Degree 
d 


Branching pattern 
above singular points 


Covering 
composition 


hypergcom. 


Hcun's 


(1/2, a, /?) 


(1/2, a, 2a, 3/3) 


3 


2+1=2+1=3 


indecomposable 


(1/2, a, f3) 


[a, a, 2a, 4/3) 


4 


2+1+1=2+2=4 


2x2 


(1/2, a, /?) 


(a, 3a, /3, 3/3) 


4 


3+1=2+2=3+1 


indecomposable 


(1/2, a, (3) 


(2a, 2a, (3, 3/3) 


4 


2+2=2+2=3+1 


no covering 


(1/2, a, /?) 


(2a, 2a, 2(3, 2(3) 


4 


2+2=2+2=2+2 


various 2x2 


(1/3, a, (3) 


(a, 2a, /3, 2(3) 


3 


2+1=3=2+1 


indecomposable 


(1/3, a, (3) 


(a, a, a, 3/3) 


3 


1+1+1=3=3 


indecomposable 



Table 1. Hypergeometric-to-Heun transformations with 2 continuous parameters 



transformed equation. We have at least c? + 2 singular points by Lemma 12.11 and we wish 
the transformed equation to be the Heun equation with four singularities, hence, d + 2 < 4 
leading to d < 2. If the degree of the covering is d = 1, we have the fractional-linear 
transformations which do not change the number of singular points. If d = 2, we have a 
quadratic transformation (a, (3, 7) < — (2a, (3, (3, 27). These transformations are discussed in 
Subsection 12.41 in full generality. 

We obtain new transformations with two continuous parameters if some points above 
{0, 1, 00} C are non-singular. (Otherwise we would have a specialization of a pull-back 
transformation of the general hypergeometric equation.) Therefore we assume d > 2 and 
restrict a local exponent difference of Hi to the value 1/k, where A; is a positive integer. 
We also assume fc > 1 because of the remark following Lemma 12.31 with a reference to 
|12[ Lemma 5.1]. The transformed equation has at most [d/fcj non-singular points above 
{0, 1, 00} C Pi, therefore it has at least d + 2 — [d/k\ singular points by Lemmas 12.11 and 
12.21 We want to have d+2 — [d/fcj <4. The left-hand side can only decrease when we drop 
the floor rounding. After ignoring the floor rounding and dividing the inequality by d we get 
the Diophantine inequality 

We have the following solutions with fc > 1, d > 2: 

(fc,d) e {(2,3), (2,4), (3,3)}. 

For each possibility (d, fc) we can write down possible branching patterns for ip such that 
the supposed transformed equation would have at most 4 singular points. The possible local 
exponent differences, the degree of the covering and branching patterns are presented in the 
first four columns of Table [1] 

Next we have to determine possible rational functions ip{x) with the branching patterns 
in the fourth column of Table [H In general, a covering with a given branching pattern may 
not exist, or there might be several different coverings with the same branching pattern. All 
our particular coverings appear in [12l Table 1] except 2 + 2 = 2 + 2 = 2 + 2, but the latter 
branching pattern determines the quartic Heun-to-Heun transformations of Subsection 12.51 
Therefore all rational functions tpi^x) are known in principle. Because the possible degree d 
is so small, one can recompute all (p{xys using a straightforward ansatz method with unde- 
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termined coefficients. In Subsection 13.41 we demonstrate a more efficient metliod introduced 
in [13l Section 3]. 

Information about existence and decomposability of tlie coverings (p{x) is given in the 
fast column of Table [H There is at most one covering for each branching pattern, up to 
Mobius transformations. 

With a possible rational covering if determined, one can choose the local exponents of Hi , 
determine the prefactor 0{x) and possible hypergeometric-to-Heun transformation formulas 
as described in Subsection 12.21 The accessory parameter of the pull-backed Heun equation 
is determined from a candidate two-term identity using Lemma 12.51 

In the following subsections we present the possible hypergeometric-to-Heun transforma- 
tion formulas up to fractional- linear transformations Q, (I6|)-(l7]). To get possible transfor- 
mations between the Heun and hypergeometric functions with one continuous parameter, we 
would restrict 2 local exponent differences of Hi to 1/fc, l/£ with k, £ positive integers, and 
would similarly get the diophantine inequalities d— [d/fcj — [d/£\ < 2 and 2/d+l/k+l/£ > 1. 
In particular, the degree bound for transformations of Hi with the local exponent differences 
(1/2, 1/3, a) is d < 12. We leave the study of these transformations for the next paper. 



3.1. The transformations (1/2, a, (3) < — (l/2j ct,2a, 3/3). Up to fractional-linear transfor- 
mations, we have the following identities: 




— 2F1 



x{4:x - 3)' 



Ax 
1 - — 



2F1 



x{Ax - 3)^ 



= 2F1 



— 2F1 



1 



a 



3x 



a, b 

-2a 



x{Ax - 3)' 
x'^{x + 3) 



— 1 2F1 



a, b 
a + b + 



(4 - 3x)2 



(23) 
(24) 

(25) 
(26) 
(27) 



where qi = {9ab + 3a + 3& — l)/4, q2 = 6a^ + 6ab — a. The five formulas represent the five 
non-equivalent choices for the local exponent difference at x = 0. The choices for the local 
exponent at a; = are, respectively: 1/2, —1/2, a, 2a, 3/3. The first two identities are related 
by Lemma [231 The arguments of the first four transformations are polynomials. Notice that 
the cubic argument in ([25]) is the same as in ([23|) -([24 |) even if the fiber for .t = is different. 
However, the branching pattern in both fibers and the branching order for a; = is the same, 
so the same configuration of the singular points a: = 0, a; = l,a; = cxDis possible (even if 
the local exponents at the respective points are different). The argument in ([26]) is related 
to a;(4a; — 3)^ by the affine transformation a; 1-^ (x -I- 3)/4, giving us other point as a; = on 
essentially the same covering. 
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3.2. The transformations (1/2, a, (3) < — {a, a, 2a, 4/3). Up to fractional-linear transfor- 
mations, we have the following identities: 



Hn 



1 

2 

8ab 



Hn 



Hn 



-1 


-1 




2a 



46 



4a, 46 
f i;a + 6- 
4a, 46 

- 26; a + 6 - 

4a, 26 

- 1:2a -6 



1 



— 2F1 
= 2F1 



a, 6 
a + b+^ 

a, 6 
a + 6+ i 

-2a 



- 16x{x~ l){2x- if 

4x^(1 - .T^ 



2F1 



a, a + 7 
6+i^ 



(x2 - 2)2 



(28) 
(29) 

(30) 



There are indeed three non-equivalent choices for the local exponent difference at x = 0, 
which are a, 2a, 3/3, respectively. The pull-back covering is a composition of two quadratic 
transformations: (^^,a,P) {a,a,2P) (a, a, 2q;, 4/3). Accordingly, the three identities 
are compositions of, respectively, (flS]) and pS]). (flB]) and ([T6|) . or (fTS]) and (flB]) . 

3.3. r/ie transformations (1/2, a, /3) <--^ (a, 3a,/?, 3/3). Up to fractional-linear transfor- 
mations, we have the following identities: 



Hnf " 



Hn 



8 

/I 



\q2 



4a, 3a -|- 6 
3a -I- 36 — ^; a + 6 

4a, 3a -|- 6 
a + 6 + i; 3a + 36 



= 1- 



= (1 



8x 
8x) 



2F1 



— 3a 17 



a, 6 
a + b + 
a, 6 



64^3(1 _ x) 

(9 - 8a;)3 
64a;(l - x)^ 
(l + 8x)3 



(31) 
(32) 



where qi — 9a^ + 9a6 — 3a/2, (72 = 3a^ — 5a6 + 3a/2. The choice between a and /3 for the local 
exponent difference at a; = gives identities related by fractional-linear transformations, just 
as the choice between 3a and 3/3. Hence we have only two transformation formulas. 

3.4. The possibility (1/2, a, (3) {2a, 2a, (3, 3/3). We demonstrate here that there is no 
covering for the branching pattern 2-t-2 = 2-f2 = 3-|-l. The three fibers are assumed to 
be z = 0, 1, 00, of course. We choose the points above 2; = 00 to be a; = and x — 00. 
(Normally, we would put the point x ~ above z = 0, so that we would have the rational 
function for transformation formula immediately. But for a negative result we have no reason 



to do that.) Suppose that the two points above z — are the roots of x^ 



b, with a, 6 



to be determined. Then the rational function defining the covering should have the form 

jx' + ax + bf 
ip{x) = C\ , 

X 

for some constant Ci . The branching pattern above z — I implies that there are two points 
above z = 1. To compute the branching points we need to compute (p'{x) = 0. We have 

, {x'^ + ax + b){3x'^ + ax - b) 

(fi (X) — Gi . 

x"^ 

The points above z = 1 must be the roots of 3x^ 4 
to a constant multiple, to 

1 - ipix) = C2 ^ 



ax — b. We must have 1 — ip{x) equal, up 
\-ax - 6)2 

X 
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The derivative of this expression must vanish at the same branching points above z = and 
2 = 1 as well. The derivative is equal to 



(3x2 



ax 



b) 



We have to conclude that the polynomials 



-b and 9x^ +ax + b must be proportional, 



which is possible only in the degenerate case a = b = 0. Hence a contradiction. 



3.5. The transformations (1/2, a, f3) 
formations, we have one identity: 



Hn 



-1 




4a, 2a -26+1 
2a + 26: 2a -26+1 



(2a, 2a, 2/3, 2/3). Up to fractional-hnear trans- 

4x2 



(l-x2) 



2\ — 2a TTi 
2ri 



a, b 
a + 6 + 



(x2 - ly 



(33) 



as the choice of the local exponent differences 2q or 2/3 for x = gives equivalent identities. 
The identity is a composition of (fH|) and In fact, the pull-back transformation can be 
expressed in several ways as a composition of quadratic transformations: 

(1/2, a, (3) ^ (a, a, 2(3) ^ (2a, 2a, 2/3, 2/3), 

(1/2, a, /3) ^ (2a, /3, /3) ^ (2a, 2a, 2/3, 2/3), 

(1/2, a, /3) ^ (1/2, 1/2, 2a, 2/3) ^ (2a, 2a, 2/3, 2/3). 

In the first two cases the second arrow denotes a quadratic Heun-to-hypergeometric trans- 
formation. In the last case, we have a quadratic Heun-to-Heun transformation. The first two 
compositions imply a hypergeometric formula like 



a, 6 
2a 



x(2-x) (1-x) 



-& / 2a - 6,6 

2-ri 

a 



(34) 



4(x- 1) 

which is a quadratic-quadratic transformation with two free parameters. This kind of trans- 
formation is presented in [21 p. 128-130]. 

3.6. The transformations (1/3, a, /3) < — (a, 2a, (3, 2/3). Up to fractional-linear transfor- 
mations, we have the following identities: 



Hn 



Hn 



9 

/ 8 



9 

V92 



3a, 2a + 6 
a + 6+ i;2a-26 

3a, 2a + 6 
2a + 26- i;a + 6- 



(1 



1 



x)-2%Fi 



9x 



2F1 



a, 6 
a + b+l 

a, 6 
a + 64 



x(x-9)2 
27(x - 1)2 

27x2(x- 1) 
(8 - 9x)2 



(35) 



(36) 



where qi — 18a'^ — 9ab + 6a, q2 — 4a2 + 4a6 — 2a/3. The choice between a and /3 for the local 
exponent difference at x = gives identities related by fractional-linear transformations, 
just as the choice between 2a and 2/3. Hence we have only two transformation formulas. 

3.7. The transformations (1/3, a, /3) < — (a, a, a, 3(3). Up to fractional-linear transfor- 
mations, we have the following identities: 



Hn 



Hn 



—UJ 

3(1 — a;)a6 

UJ+1 

3{lu + 2)a6 



3a, 36 
'6+i;a+6+i 

1 

3 
2 
' 3 



3a, a + 6 - 
36: 2a - 6 



= 2F1 
= (1- 



a, 6 
a + 64 



3(2w + l)x(x-l)(x + w) , (37) 



a, a + 
6+1 



{x-Lu-2)^. 



(38) 
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where u) is the root oi lo'^ + uj + \ ~ {). The choices for the local exponent at a; = are 
a and 3/3. To relate the argument in (j37p to [71 formula (3.6a)], note that 

3(2w + l)x{x ~ l){x + cj) = 1 - (1 - + 2)xf. 

The other two-parameter coverings of [7] are represented in Subsections 13.11 and 13.21 above. 
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